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Abstract There is a general consensus about the fact that the magnetar scenario pro- 
vides a convincing explanation for several of the observed properties of the Anoma- 
lous X-ray Pulsars and the Soft Gamma Repeaters. However, the origin of the emis- 
sion observed at low energies is still an open issue. We present a quantitative model 
for the emission in the optical/infrared band produced by curvature radiation from 
magnetospheric charges, and compare results with current magnetars observations. 



1 Introduction 

The Anomalous X-ray Pulsars and the Soft Gamma Repeaters (AXPs and SGRs) 
are peculiar X-ray pulsators which exhibit erratic phases of bursting/outbursting 
activity (e.g. Mereghetti 2008). There is now a wide consensus that these sources 
indeed host an ultra-magnetized NS, or magnetar (Duncan & Thompson 1992), with 
a magnetic field B w 10 14 — 10 15 G, well in excess of the quantum critical field Bq ~ 
4.4 x 10 13 G. 

The magnetar scenario provides a quite coherent interpretation of the observa- 
tional properties of SGRs/AXPs, in particular of both their bursting and persistent 
emission at X/y-ray energies (Wood & Thompson 2006). In particular, it has been 
suggested that the magnetosphere of a magnetar is twisted (VxB^ 0), in which 
case the density of magnetospheric currents is substantial and the observed soft X- 
ray spectrum (< 10 keV) can be explained in terms of resonant cyclotron scattering 
(RCS) of thermal surface photons onto mildly relativistic particles flowing along 



Silvia Zane 

Milliard Space Science Laboratory, University College of London, Holmbury St Mary, Dorking, 
Surrey, RH5 6NT, UK, e-mail: sz@mss l . ucl . ac . uk| 

Luciano Nobili, Roberto Turolla 

Dept. of Physics, University of Padova, Via Marzolo 8, 35131 Padova, Italy, e-mail: 
|luciano . nobili @pd . inf n . it , roberto . turolla@pd . inf n . it"| 

1 



2 



The magnetar emission in the IR band 



the closed field lines (see Zane et al. 2009 for a recent application, and references 
therein). 

Observations at IR/optical wavelengths led to the discovery of faint (K ~ 19- 
21), and in some cases variable IR counterparts to several AXPs/SGRs (see again 
Mereghetti 2008). In two sources optical pulsations at the X-ray period were found 
(Kern & Martin 2002, Dhillon et al. 2009). However, the origin of the IR/optical 
emission is not clear as yet. The proposed scenarios involve emission/reprocessing 
by a fossil disk (see e.g. Perna et al. 2000), or non-thermal emission from the in- 
ner magnetosphere (Eichler et al. 2002, Beloborodov & Thompson 2007, hereafter 
BT07). Although this latter option appears promising, no detailed computations 
have been presented as yet. In this paper we reassess this issue, and provide a quan- 
titative estimate of the IR emission expected from curvature radiation. 



2 Pair production in the inner magnetosphere 

In the inner magnetosphere, the magnetic field is so strong that the conversion 
of high-energy photons produced by RCS into pairs may give rise to a quasi- 
equilibrium configuration in which the e ± Lorentz factor is locked at about the 
threshold value required to ignite the pair cascade (50 < y < 1000, BT07). In order 
to estimate the spatial extension of this region, we proceed as follows. We assume 
that primary photons have a typical energy £ ~ 1 keV, and the magnetic field scales 
as B ~ Bp(R/Rns)~ 3 , where B p ~ \QB q is the polar field and R^s is the star radius. 

We assume that the magnetosphere is force free. Charges, flowing along B, are 
accelerated by the parallel electric field which is provided by the decay of the twist, 
d(B 2 /87r)/dt = —E\\j. A thermal photon scatters where its energy in the parti- 
cle frame, 7(1 — /3cos0)e, equals the local cyclotron energy, = m e c 2 {B/Bo) 
(here j3 = v/c ~ 1 and 9 is the angle between the incident photon and the mag- 
netic field), i.e. when the particle Lorentz factor is y = y res ~ (m e c 2 /e)(B/Bq). 
Because the potential drop at the endpoints of a flux tube (of length L) is enormous, 
e<Po > 10 6 ra e c 2 (see below), the electron Lorentz factor, which is ~ 1 at the star 
surface, monotonically increases up to a value which is much higher than y res ,max ~ 



Fig. 1 Geometry of the cur- 
vature emission. See text for 
details. 
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500(1 keV/e) (B p /Bq) ps 5000. This ensures that the resonant condition (Yres = y) 
is met after the electron travelled a typical distance X acc ,res = Y\dx/d(y^)\ ~ 
YresL/ (e<t>o/m e c 2 ) <C L. On the other hand, this is not a sufficient condition for RCS 
to occur: even if they can be accelerated up to y res , electrons are able to scatter only 
if X < X c ,ccj- es , where X - (e/m e c 2 )/ (n ph a e ff) ps 10~ 3 (r /IkeVy 2 (R/R NS ) 2 cm is 
the particle mean free path for resonant scattering (here is the density of thermal 
photons and a e /f = 3kGt/40Cf; e.g. Dermer 1990). This translates into: 

m e c 2 \B Q J\R NS J \lkeVj Vl0 6 cm y ' 

We approximate the left hand quantity with the potential drop which produces the 
conduction current in a ID, relativistic double layer (where j oc <J>~; Carlqvist 1982): 

£^ 2xl0 » 1 /^(-Ly / VAy"V4_v (2) 

m e c 2 V m e \B Q J \R NS J V lo6cm / 

where m + is the mass of the positive charge. If the magnetosphere is populated by 
e~ /e + pairs then m + = m e , and for a polar magnetic field <~ IOBq the inequality 
given by Eq. <£TJ is satisfied up to B/Bq > 0.05, or R/Rns < 6. 

The upscattered photon produced in the region B/Bq > 0.05 has an energy 
e ' ~ Tre.s e /(1 + Jres£/m e c 2 ) and may convert into a e ± pair in the strong magnetic 
field provided that the threshold condition e' > 2m e c 2 / sin 0' is met, where 9' is the 
angle of the scattered photon with the magnetic field. For e <~ 1 keV, pair creation 
requires y res > 30 (i.e. B > Q.QSBq). Since electrons are relativistic, the scattered 
photon initially moves parallel to B and a pair can be produced only after a large 
enough pitch angle has built up. Moreover, pair production is efficient only if its 
characteristic length-scale is less than the dimension of the circuit, i.e. a± > 1 /Rns 
where <x± is the absorption coefficient for ordinary and extraordinary photons (see 
Baring 2007). Accordingly, the limiting value of magnetic field above which the 
process is effective decreases with increasing photon energy, and it is B <~ 0.05Bg 
for pairs created near threshold. 

In summary, in the region B > Q.QSBq a quasi-equilibrium configuration is 
reached with a pair multiplicity ~ L/X aC c.res of a few. Screening of the electric field 
limits the potential drop to e<Po/meC 2 ~ Yres ~ 5QQ(B/Bq) and the maximum e ± 
Lorentz factor is Yres- Charges undergo only few scatterings with thermal photons, 
but they loose most of their kinetic energy in each collision. A steady situation 
is maintained against such Compton losses because e~/e + are re-accelerated by 
the electric field before they can scatter again. RCS may occur also at larger radii 
provided that the charges Lorentz factor is limited to ~ y rev , which decreases as 
B ~ The pairs energy depends on the electrostatic potential (in turn fixed by 
conduction current), and also on the efficiency of the radiative (Compton) drag. If 
pairs with y ~ y res are injected in the external region (B/Bq < 0.05), the circuit is 
likely to behave much differently from a double layer, allowing the current to be 
conducted with only a small potential drop (see also Beloborodov, this volume). 
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3 IR/Optical emission 

In order to estimate the amount of curvature emission from the inner magnetosphere, 
as seen by a distant observer, we use a simple geometrical model. We approximate 
the twisted field with a dipole, B = B p (R/Rhs)~ 3 (cos 9, sin 0/2,0), while the cur- 
rent is still taken to be j = (c/4n)\W x B\ tW i St « B(R/Rns)^ 1 ■ The particle density 
is then n± = B(R, 9) (R/Rns)^ 1 / (4fte), for a pair velocity ~ c. With reference to 
a generic flux tube labelled by 9q (0 < 9$ < n/2), curvature photons will reach 
the observer if they are emitted in the volume A V = XR t sin 9 \, s \ dZ/d9o\A9o/y 2 , 
around the point P of cartesian co-ordinates (X, Z) where B is parallel to the LOS 
(FigJTJi. Here 9 t, s is the viewing angle,i.e. the angle between the line-of-sight (LOS) 
and the magnetic axis. In general there are two emitting regions along each flux tube 
one above the other below the magnetic equator, because of electrons and positrons, 
and for any 9q there are two azimuthally symmetric flux tubes. 

The (polarization averaged) CR spectral power per particle can be expressed as 

P (e) = {E/ec) (e/e c )f" /ec K 5/3 (x)dx / [^xdx^K 5 / 3 {x!)dx!], where K 5/3 (x) is 

the modified Bessel function, e c = (3/2)/zcy 3 /R c , E = (2/3)e 2 cf /R~, and R c is 
the curvature radius of the field line. The monochromatic luminosity emitted by the 
volume AV is AL = p(e)n±AV, and, since radiation is collimated in the solid angle 
AO. = n/y 2 , the observer at infinity measures a luminosity 4-nAL/ AO.. The total 
received power is then obtained summing the individual contributions over all the 
flux tubes, assuming the pair Lorentz factor is locked at y res . 

Fig.|2](left panel) shows the computed CR spectrum for B p = 4.3 x 10 14 and dif- 
ferent viewing angles 9 (> i, s , together with the optical/IR luminosities observed for a 
sample of AXPs. It should be noted that all these sources typically show variability 




Fig. 2 Model spectra for different values of 9„ bs and B p = 4.3 x 10 14 G. The XMM-Newton X- 
ray spectrum of 1RXS J 1708-4009 is from Rea et al. (2008; red solid line). The AXPs IR/optical 
data are from Duncan & van Kerkwijk 2005 (4U 0142+614 and IE 1048-5937) and Mignani et 
al. 2007 (XTEJ1810-197 and 1E2259+586). The adopted distances for de-reddening are 5 kpc 
(4U 0142+614, 1RXS J1708-4009), 3 kpc (IE 1048-5937, IE 2259+586), 4 kpc (XTE J1810-197), 
8.5 kpc (IE 1841-045). Left: Spectra computed for incoherent curvature emission. Right: Same as 
in the left panel, but accounting for particle bunching. 
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in all energy bands. The data shown in Fig. [2] were not selected according to any 
particular criterion and are just meant to be representative of the IR emission of 
magnetar candidates. As it appears from Fig. [2] model spectra severely underpre- 
dicts optical/IR data. Basically, this is due to the fact that the particle density is too 
low for incoherent CR emission to produce the required power. The point has been 
already noted, on the basis of a qualitative analysis, by BT07 who concluded that 
particle bunching would be likely present, providing the required amplification of 
CR emission. 

Many models have been suggested to explain the origin of the interactions which 
push particles together and can lead to the formation of bunches of charged particles 
localized in phase-space, but the most promising explanation seems to be connected 
with plasma instabilities, like the two-stream (electron-positron/electron-ion) insta- 
bility. If N particles are contained in a bunch of spatial scale l B , they radiate like a 
single particle of charge Q = Ne. This collective emission produces an amplifica- 
tion of the radiated power by a factor N with respect to the standard, non-coherent 
emissivity by N individual charges. The details of the emission depend on the shape 
of the bunch. In the simplest, one-dimensional case the power per unit frequency 
emitted by a bunch is Pb(v) ~ N 2 p(v) sin 2 (tzv/vi)/(7Zv/vi) 2 = N 2 3 g (v) (Sag- 
gion 1975). The main source of uncertainties lies in the size of the bunch, and 
here we assume that Zg is associated to the local plasma frequency, lg ~ c/v pe = 
%m e c 2 /{je 2 n e ) (e.g. Lesch 1998). Since a single bunch emits a spectral power 
(n e lg) 2 3? (v) and there are Nb ~ 4V//| bunches in each emitting volume, the total 
spectral power radiated within AV is ALb = NB{n e ll) 2 £P(v) = (n e l B )AL. We stress 
that it has to be N = n e \\ > 1 for the process to be effective, meaning that only 

1/3 

emission at frequencies v < V/ < V co = n e c is efficiently amplified. Where N <\ 
emission is from individual particles. Finally, because of the strong absorption be- 
low the plasma frequency, a low-energy cut-off is present around v pe . 

We have recalculated the emergent spectrum including particle bunching (Fig. [2] 
right panel), finding that coherent curvature radiation can indeed produce sufficient 
IR/optical emission to account for the observed one. Given the limitations of our 
model, we are not in a position to attempt any fit of the data. What we aim to, 
instead, is verifying that magnetospheric CR emission is able of reproduce the gross 
characteristics of the observed low-energy emission, its energetics in particular. The 
predicted spectrum vary with the viewing angle, and the luminosity is higher when 
the star is viewed nearly pole-on. The curvature spectrum is far below the observed 
one in the soft X-rays, indicating that a different mechanism is necessary to account 
for the X-ray emission below ^10 keV. In the twisted magnetosphere scenario this 
is RCS onto mildly relativistic pairs populating the external magnetosphere. 



4 Discussion and Conclusion 

Despite the magnetar scenario is now regarded as the most likely interpretation for 
the observed properties of SGRs and AXPs, and the many investigations aimed at 
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producing detailed models for the spectral and timing properties of magnetars, a 
number of key issues are still unresolved. In particular, the origin of the low-energy 
emission (in the IR/optical band) of a magnetar is still under debate. In this paper we 
gave a quantitative estimate of the IR/optical emission, under the assumption that it 
arises from curvature radiation from pairs in the inner magnetosphere. A comparison 
with observations of AXPs shows that the predictions of the model are in general 
agreement with observations, provided that curvature radiation is not coherent and 
a bunching mechanism is at work at long wavelengths. 

Our model is based on a number of simplifying assumptions. For instance, the 
computation is based on a globally twisted dipolar magnetosphere while in some 
magnetars, the transient AXPs in particular, the twist might affect only a limited 
bundle of field lines close to a magnetic pole (Beloborodov 2009) . If the inner 
part of the magnetosphere is (or becomes) untwisted, no currents are present that 
can produce the IR/optical curvature radiation. Also, we used a simple double layer 
model (eq. 1) for the linear accelerator, which may not be valid in presence of intense 
pair creation (BT07). Furthermore, and more important, the details of the charges 
motion in the external region of a twisted magnetosphere, which are essential ingre- 
dients in providing a model for the high energy emission in the soft and hard X-ray 
band (~ 0.5-200 keV), are still unclear. In order to reach firmer conclusions about 
the entire multi-wavelength spectrum a more detailed study of the magnetosphere is 
required, and will be matter of future work (see also Beloborodov, this volume). 
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